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Abstract. We address the problem of the foundation of generalized ensembles in 
statistical physics. The approach is based on Boltzmann's concept of orthodes. These 
are the statistical ensembles that satisfy the heat theorem, according to which the heat 
exchanged divided by the temperature is an exact differential. This approach can be 
seen as a mechanical approach alternative to the well established information-theoretic 
one based on the maximization of generalized information entropy. Our starting point 
are the Tsallis ensembles which have been previously proved to be orthodes, and have 
been proved to interpolate between canonical and microcanonical ensembles. Here we 
shall see that the Tsallis ensembles belong to a wider class of orthodes that include the 
most diverse types of ensembles. All such ensembles admit both a microcanonical-like 
parametrization (via the energy), and a canonical- like one (via the parameter (3). For 
this reason we name them "dual" . One central result used to build the theory is a 
generalized equipartition theorem. The theory is illustrated with a few examples and 
the equivalence of all the dual orthodes is discussed. 
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1. Introduction 

Around 20 years ago a generalization of the standard Gibbs-Boltzmann statistics was 
proposed by C. Tsallis The generalization scheme proposed by Tsallis was based on 
the adoption of a generalized Shannon's informational entropy, which, when maximized 
under the constraints of normalization and average energy, leads to power-law statistics. 
Since then, the corresponding generalized thermostatistics, also named non-extensive 
thermodynamics, has proved to be a very powerful tool of investigation within the most 
diverse fields of physics. 

In two previous works (2j [3] we have proposed an alternative mechanical approach 
to study the foundation of the so-called Tsallis ensemble. Those studies revealed 
two interesting facts. First, Tsallis ensembles are exact orthodes [2]. According to 
Boltzmann's original definition (see [lj for a modern exposition) these are the ensembles 
that satisfy the "heat theorem" 

= exact dif ferential (1) 

where 5Q = dU + PdV and T are calculated from the ensemble average of properly 
chosen functions of the microscopic state (i.e. the 6N-dimensional phase space 
vector). In this sense it is said that the orthodes provide good mechanical models 
of thermodynamics |4j. Second, not only the canonical ensemble is a special case 
of Tsallis ensemble, but the microcanonical is also. That is, the family of Tsallis 
ensembles parameterized by the extensivity index q interpolates between canonical and 
microcanonical ensemble [3J. 

We now have a continuous family of ensembles, namely the Tsallis ensembles, 
parameterized by the extensivity index q that generalizes both the canonical and 
microcanonical ensembles. All these ensembles are orthodes. Here we shall generalize 
further and show that said family is a subset of a more general class of orthodes which 
we shall call dual orthodes. The main ingredient that will be used to develop the 
generalization is the duality property which characterizes the Tsallis ensemble. By 
duality we mean that the ensemble can be parameterized either through the average 
energy (microcanonical-like parametrization) or through the parameter {3 (canonical- 
like parametrization) . In section [2] we provide the definition of orthode and we review 
the content of Refs. [2, 3J. We move to propose the generalization in Sec. [3J where 
the concept of dual statistics is introduced and a generalized equipartition theorem is 
presented. Section H] shows that the standard ideal gas thermodynamics follows from 
any dual statistics. In section [5] we illustrate the concept of dual ensembles by studying 
a few particular cases, which include gaussian and Fermi-like ensembles. Section [6] gives 
the concluding remarks. 
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2. The Tsallis orthodes 

This section is devoted to review the results of Refs. [21 E]. These concern the foundation 
of Tsallis statistics from the viewpoint of orthodicity. The approach can be considered 
a mechanical approach, as opposed to the standard information theoretic one used by 
Tsallis pp. The concept of orthode was first proposed by Boltzmann, but did not meet 
the same success of other ideas developed by him, such as the counting method or the 
H-theorem. Nonetheless the method provides of a very straightforward way to asses 
whether a given statistical ensemble provides a mechanical model of thermodynamics. 
This method was used by Boltzmann to give a theoretical basis for the canonical and 
microcanonical ensemble in statistical mechanics. Now we know that Tsallis ensembles 
also provide good models of thermodynamics as well [2]. The idea on which orthodes are 
based is quite simple. Consider a generic statistical ensemble of distributions, p(z; Ao, Aj) 
parameterized by one "internal" parameter Ao and a given number of "external" 
parameters Aj i > 0. The symbol z denotes the phase space point vector. We assume the 
system to be Hamiltonian with z = (pi,j9 2 , ■■■,p n ,Qi,Q2, ■■■,Qn)- The Hamilton function 
is assumed to be of the form H = K(p) + <£>(q; Aj), i > 0. The Hamiltonian depends 
explicitly on the external parameters, which are also called "generalized displacements" . 
For sake of simplicity, we consider only one external parameter Ai = V. This may be, for 
example, the volume of a vessel containing the system. The kinetic energy is assumed 

2 

to be of the form K(p) = j-. Then define the macroscopic state of the system by the 
set of following quantities: 

U p =(H) p "energy" 

T p = "doubled kinetic energy per degree of freedom" , 

V p = (V) "generalized displacement" 

P p = ( — If?) "generalized conjugated force", 

Now let the parameters Ao, V change of infinitesimal amounts, and calculate the 
corresponding change in the macroscopic state. If the state changes in such a way 
that the fundamental equation of thermodynamics, namely 



dU p + P p dV p 
T 



exact dif ferential, (3) 



holds, then we can say that the ensemble provides a good mechanical model of 
thermodynamics. Equation ([3]) is also known as the heat theorem and expresses the 
fact that the heat differential SQ admits an integrating factor (T _1 ), whose physical 
interpretation is that of reciprocal temperature (namely doubled average kinetic energy 
per degree of freedom). Correspondingly the function which generates such exact 
differential is interpreted as the physical entropy of the system. The ensembles that 
satisfy Eq. ([3]) were called orthodes by Boltzmann [5] (see also [Ij and [6] for modern 
expositions). Boltzmann proved that microcanonical and canonical ensembles are 
orthodes, and placed this fact at the very heart of statistical mechanics. For example, 



Thermodynamics with generalized ensembles: The class of dual orthodes 4 

the canonical ensemble is parameterized by one internal parameter, usually indicated 
by the greek letter (3, and the external parameter V: 

-f3H(x;V) 

»^W = lWV) (4) 

where 

Z(P,V) = J dze- pH W\ (5) 

The following function: 

S e (/3, V) = (3U C {(3, V) + log Z(/3, V) (6) 

generates the heat differential, therefore is the entropy associated with the canonical 
orthode. To prove that, we have to calculate the partial derivatives of S c : 

where the symbol < • > c denotes average over the canonical distribution and the 
state definition ([2} has been used. Further according to the canonical equipartition 
theorem 

Therefore, by comparison with ^2} T c = 4. Combining all together we get: 

^f dU ^a^9U \ dU c + P c dV c nm 

dS = p [w dp + W dv + Pcdv ) = — t c • (10) 

The microcanonical ensemble is parameterized by (U, V): 

^.vj-ffi^Qi (U) 

where 5 denotes Dirac's delta function and 

Q(U,V) = J dz5{U-H{z;V)) (12) 

denotes the density of states or structure function |?]. As shown in Ref. [6] the 
corresponding entropy is given by 

S mc (U, V) = log J dz6(U - H(z; V)) (13) 

where 9 denotes Heaviside step function. In this case the proof is based on the 
microcanonical equipartition theorem [7]: 

/ BH\ = W I V 1 

mc v></L mv) ( > 
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where 



where < • > mc denotes average over the microcanonical distribution fill I) and 

$(U,V) = Jdz6{U-H{z; V)). (15) 

A proof that the heat theorem is satisfied by the microcanonical ensemble appeared in 
i- 

The Tsallis ensembles of indices q < 1 are dual ensembles (2j |3] , namely they can 
be parameterized either through (/3, V) or (U, V): 

p q (z;U,V) = , (16) 

iV <? 

N g (D, V) = Jdz[l-(1- q)/3(H(z- V) - U)]^ . (17) 

The symbol "□" has to be replaced by either U or (3 according to the parametrization 
adopted. In the (U, V) parametrization, one first fixes U and then adjusts (3 in such a 
way that: 

U = (H) pq , (18) 

In this parametrization Eq. (JTHJ) defines the function /3(U,V). In the ((3, V) 
parametrization, f3 is fixed instead and U is adjusted accordingly, in such a way that Eq. 
( [18]) defines the function U(/3, V). As shown in Ref. [2], no matter the parametrization 
adopted the ensembles (fl6l) are orthodes and the corresponding entropies are given by: 



S [,1 (D, V) = log J dz [1 - (1 - q)[3{H{z; V) - U)}^ (19) 
The proof of orthodicity is based on the corresponding Tsallis equipartition theorem: 

m 



dH\ lN q (D,V) 



where: 

JV q (D, V) = J dz [1 - (1 - q)P(H(z; V) - U)]^ (21) 

The proof of orthodicity of Tsallis ensemble appeared in Ref. [2] . 

The fact that Tsallis ensembles are orthodes reveals that they are as well founded as 
the microcanonical and canonical ensembles, at least as far as providing good mechanical 
models of thermodynamics. The work of Ref. [3] has shown that the family of 
Tsallis ensembles parameterized by the index q interpolates between canonical and 
microcanonical ensembles. In fact, the former is recovered in the limit q — > 1, and 
the latter is recovered in the limit q — > — oo. This interpolation is not limited to the 
distributions only but also to the corresponding entropies and equipartition theorems. 
Physically the Tsallis statistics describes the situation of a system in contact with a 
finite heat bath, namely, a bath with a given finite heat capacity [U [9]: 

C v = -r^— (22) 
1 - q 
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As q goes to 1, the heat capacity goes to infinity and correspondingly one obtains the 
canonical ensemble (the system is thermalised) . As q goes to — oo, the heat capacity 
goes to zero and one recovers the microcanonical situation (the system is thermally 
isolated) [3]. 



3. Generalization: Dual orthodes 



So far we have seen how the Tsallis ensembles of index q are orthodes with a special 
duality property We have also mentioned that the microcanonical and canonical 
ensembles are two special instances of Tsallis ensembles. Now we shall see that the 
property of orthodicity can be proved for a general class of ensembles which share with 
the Tsallis ensemble the fact that both U and (5 appear explicitly in their expression. 
They can be considered as parameterized by either U, V or (3, V, depending on which 
parameter is kept fixed and which one is adjusted, in such a way that U =< H >. We 
shall call these ensembles dual ensembles. 

The generalization is purely formal, in the sense that we shall assume that all 
integrals and derivatives written exist. Further we shall assume that, for given U, V (or 
P, V), the equation U =< H > admits a solution P(U, V) (or U(f3, V)) which is of class 
C 1 . These are conditions that must be checked, a posteriori, on a case by case basis, 
depending on the explicit form of the Hamiltonian and of the distribution. Thus, let us 
consider a generic ensemble of the form: 

Where G and j3 are assumed to be differentiable functions of U, V. This means that we 
are adopting the U, V parametrization, alternatively we could have adopted the /3, V 
parametrization. The values of G and j3 are fixed by the two constraints 

G(U, V) = J dzf[(3(U - H(z; V))} (24) 

and 

U = J dzH{z;V)p{z;U,V) (25) 

Throughout this section the symbol < • > p denotes average over the dual distribution 
( I2"3"j) . Consider now a differentiable function F(x) such that F(x) > 0,F'(x) = f(x), 
and define the following function: 

S p (U,V)=\o g g(U,V) (26) 

where 

Q(U, V) = J dzF[P(U - H(z; V))} (27) 

and all the integrals are extended to the definition domain T> of the distribution p. We 
shall also assume that a cut-off condition exists such that F is null on the boundary of 
the domain: 

F[(3(U-H(z;V))] Zedv = 0. (28) 
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Let us define the macroscopic state: 
U p = (H) p = U 



2<*>p 



v p = v 



(29) 



Before proving that the ensembles of the form (123]) are orthodes let us state the 
following generalized equipartition theorem: 

Theorem 1 In the U,V parametrization the average ofpi 911 

dH\ 1 G(U,V) 



is: 



Pi 



d Pi / p p(u,v)G(u,vy 



In the (3,V parametrization the average of Pi^f- is: 



P 



dH\ l Gj^V) 
d Pi / PG{J3,V) 



(30) 



(31) 



The proof is provided in Appendix A It involves writing the integral expression of the 
quantity (vijjfj , integrating by parts over p^ and using the cut-off condition fl28l) . The 

proof structure is the same as the structure of the proof of Tsallis equipartition theorem 
ofRef. 0. 

Comparing with the macroscopic state definition ( 1291) . the generalized equipartition 
theorem can be reexpressed as the following compact formula : 

ig 

(3G 

where it is intended that Eq. (130]) is used in the U, V parametrization and Eq. (f3~Tj) in 
the P, V one. 

Equation (|32l) tells that for generic dual statistics the quantity (3 might not coincide 
with the inverse physical temperature. Let us now evaluate the partial derivatives of 
the entropy function in Eq. ( 1261) : 

dS n 



(32) 



dU n 



l -jLj d .FW-H)) 

d(3 



dzf(/3(U-H)) 
Gdf3 



dU 



G 



gdu ^ H - u) » + ¥ 



T n 



(33) 



In order to obtain the last equality we used the first definition in ([291 and Eq. (|32~1) . 



dSp 
dV n 



1 d C 

-— J dzF((3(U-H)) 

dp 



1 

Q 
Gdp_ 
GdV 



dzf(P(U-H)) 



(U-H)-P 



OH 
dV 
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= ^- (34) 

In order to obtain the last equality we used the first and fourth definitions in (1291) and 
Eq. ([32]). From Eq.s (|33l) and ([34} we get: 



dE p + P p dV p 

J P ~ rp 

dE +P dV, 



= ^PJ_±^P (35) 



Therefore the differential — e -^ £ — - is exact and the entropy is given by Eq. (126|) . This 
implies that the ensembles of the form (1231) are orthodes, namely they provide good 



mechanical models of thermodynamics. In Appendix B we provide a proof that the 
heat theorem is satisfied also if the alternative j3, V parametrization is adopted. The 
proof is essentially the same as for the Tsallis case (see [2]). Thus, we have found that 
the class of orthodes, whose known representatives have been for more than one century 
only a few (canonical, microcanonical, grand-canonical and pressure ensemble jl]) is 
indeed quite vast and can include other statistics. 

3.1. Recovery of known cases 

3.1.1. Canonical The canonical ensemble is a very special case of dual orthode where 
the parameter U does not appear explicitly in the expression of the distribution. This 
case is obtained with the choice: 

f(x) = F(x) = e x . (36) 

In this case we get 

e f3(U~H) e -f3H 

P(Z; & V) = JdnePV-B) = Jd^m- (37) 



The average energy U cancels in the last term of (1371) . In this sense we refer to the 
canonical ensemble clS 8b CclSG of "hidden dual orthode". The canonical entropy is 
recovered by taking the natural logarithm of Q = J d7,e^ v ~ H ^ : 

S(/3, V)=(3U + log J dze-? 11 (38) 



Note also that, from Eq.s (Jz4j) and (J27|) . G = Q in this specific case, therefore the 
generalized equipartition theorem ( 1321) gives T p = 4. In this way the canonical 
equipartition theorem (jUJ) is recovered too. 

3.1.2. Microcanonical The microcanonical ensemble is recovered with the choice 

f(x) = 6(x) 
F(x) = 6(x) 

From the properties of the Dirac delta the distribution in Eq. ( 123]) is: 
p{z, U,V ) - j _ H)) ~ j dzS{u _ R) 
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As with the canonical case (1371) . the last term in f)39p . does not depend explicitly 
on j3, hence the microcanonical case is also a case of "hidden dual statistics". The 
microcanonical equipartition theorem ( Tl4|) is also recovered. From ( 132]) one gets: 

1 g _ J dz6(U - H) _ $ 
]3G ~ Jdz5{U-H) ~ H 

where we have used the relations 6(ax) = 6(x) (for a > 0) and S(ax) = a^ 1 5{x) 



T p = - a ^ = \ZTrr rri = 7i ( 40 ) 



3.1.3. Tsallis The Tsallis distribution is recovered with the choice 

f(x) = [l + (l-q)x]& (41) 
F(x) = [1 + (1 - q)x] T ^ 

(42) 

In this case one finds Q = M q and G = N q in Eq.s ( BID an< ^ (EZJ), so from Eq. ( 1321) the 
Tsallis equipartition theorem fTSOj) . is recovered: 

T.AiA%- (43) 



Canonical and microcanonical cases are both included in the family of Tsallis 
distributions as special cases corresponding to the values q — 1 and g = — oo [2]. 
From [21 E] we have for finite g's Af q = N q , therefore Eq. fT4"3"|) becomes: 

T M = i |g| < oo (44) 
4. Derivation of the ideal gas thermodynamics 

In the ideal g the potential energy (p(z; V) is a box potential which constrains the 

coordinates to lye in an interval of measure L = W where we assume for simplicity a 
cubic box of volume V. The box potential reduces the integration over the configuration 
space to a domain of measure V n ^ 3 , where n = 3N is the total number of degrees 
of freedom and iV is the number of particles. The Hamiltonian is purely kinetic 
H = Y^=i J 1 *- Assuming that the function /3(U, V) exists, the equation of state is 
obtained from Eq. (|34|) . namely = jr: 



o 



d n q J d n pF[P(U - H(p))\ 



8S P _ 1 d 

W ~ gov 

1 d : V n/3 J d n pF[/3(U -H(p))} 



gov 

-17 < 45 » 

from which the standard ideal gas law is easily obtained: 

P P V P = gT p (46) 

The fact that the standard ideal gas law is found to hold for any dual statistics 
should not surprise since it generalizes a result already found within the nonextensive 
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thermodynamics [TO]. Let us now focus on the form of the function (3(U, V) in the ideal 
£cas case. Using the standard change of variable [7] = e, d n p = c n ea 1 de, followed 
by the change of variable y = /3e, the condition < H >= U is expressed as: 

rv{A) 

I (A, n) = / dyy^\A - y)f(A - y) = (47) 
Jo 

where A = /3U, and y is a cut-off possibly infinite and possibly depending on A. We 
shall refer to Eq. (j!7|l as the energy constraint equation. From such equation it is easily 
inferred that, if for a given number of degrees of freedom n, a solution A n of ( )47l) exists, 
then the function (3{U, V) exists and is given by: 

P = ^ (48) 

For example, within the canonical ensemble one has (3 — In the ideal g as CRSG 
< K >= C/, so, from the definition of macroscopic state fl29|) one has 

T p = ™. (49) 
n 

Because of orthodicity we have ^ = = ~TU ^ or an y dual ensemble. This implies that 
the entropy is: 

- log U + - 
2 s 3 

Therefore the standard ideal gas thermodynamics is recovered for any dual orthode. This 
means that the canonical or microcanonical statistics are not the only necessary statistics 
that give the ideal gas thermodynamics. On the contrary ordinary thermodynamics may 
follow from non-ordinary ensembles belonging to the class of dual orthodes. 

As the generalized equipartition theorem (1521 suggests, in general the standard 
relation T p = i does not hold. For example, from Eq.s (HBl) and (j49p . one easily finds 
the following formula: 

T ' ^ ik (51) 

By comparison with the Eq. fl32l . one also deduces that, 

(52) 



S p (U, V) = - log U + - logy + const (50) 



r> n 

y 2 



G A n 

The previous formula can also be derived directly by considering the explicit expression 
of Q (we adopt the U, V representation), 

g(U,V) = c n V% dee^FlA^l + e/U)] (53) 
Jo 

where U is the cut-off energy value. Equation (1521) follows after an integration by parts 
and the definition of U from Eq. ( |29|) . 
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5. Examples 

5.1. Tsallis ensemble 

As an illustration of the theory let us first apply it to the Tsallis orthodes \2\ of indices 
q. As we will see, this is a quite special case that can be worked analytically. The 
ensembles are (we adopt the U, V representation): 



: (H(z;V)-U) 



(54) 



N(U,V) 

where for simplicity we have adopted the notation a q = y^— . In this case the energy 
constraint integrals ([17)1 can be evaluated analytically: 



/„(An) = -Q-x)^ +a «»" ,+ " /2 ' lr(Q ' +1)r ® 



r(f + « g + i) (55) 

The solutions A n>q of the equations I q (A,n) = are A nA = |, no matter the value of 
q. Therefore one finds, from Eq. (148j) . the relation T^' = which is in agreement with 
the Tsallis equipartition theorem (see Eq. (fl4|) ) 

Using Eq. (jJSJ) one can express the Tsallis ensemble, in the ideal g clS CclSG, clS 



p q (v,u,v) 



9 

1-9 



1 



(i-g)n 



iJ(z;V)' 



2 V / 



1-9 



(56) 



where the £7, 1/ parametrization has been adopted. Alternatively, adopting the /?, V 
parametrization the Tsallis ensemble would read as: 



p q (z;U,v) 



q){*-PH{*\V))Y 



Applying Eq.s 



Jdz [l + (l- g )(s-/3JT( Z ;V))] 
and (1271) gives the entropy. In the U, V representation, it reads: 



(57) 



(U, V) = 5 log V + log L n JU) + log c v 



where 



L ntq {U) = Z/f 



i+- 



dxx 2 



(1 — o)n . , 
l+ l ^ 1-x 



i 

1-9 



(5f 



(59) 



In agreement with Eq. (1501) . the dependence of the entropy on U is of the type S logt/. 
The integral L nj<? has been obtained by using the cut-off condition e = U (^1 ' 
and the change of variable x = e/U, where e denotes energy. 



n(l-q) I ' 



5.2. Gaussian ensemble 

Since the fundamental work of Khinchin [7] based on the application of the central 
limit theorem, it is known that the distribution law for a large component of a large 
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Hamiltonian isolated systems of total energy a is well approximated by a the following 
Gaussian distribution: 

(60) 



e (3{a-H) exp 



P 



2B 2 



normalization 

the quantities A\ and B 2 being defined in terms of the Laplace transforms Zi((3) of the 
structure functions Qi(x) of the system (i = 1) and the heat bath (i — 2): 

dhg Z x 



B 2 



dp 
d 2 \ogZ 2 
dp 2 



(61) 



According to Khinchin the quantity A% is a good approximation to the average energy 
U of the system ( t/ "^ / Al = where N\ is the number of particles in the system). 

Besides B 2 , the width of the distribution, can be expressed, in the case of an ideal gas 
bath (see Chapter 5, Section 22 of Ref. [7]) as: 

37V 2 



Bo 



2/3 2 



(62) 



Here N 2 denotes the number of particles in the bath. Hence the ensemble in Eq. (1601) 
can be re-expressed in the form of a dual ensemble: 



p a (z;U,V) = e 



cxp 



2(T 



G 



(63) 



where a=^fis indeed the specific heat of the heat bath, namely a plays the same role 
here as the parameter a q = y^— in the Tsallis ensembles (see Eq. (1221)). The Gaussian 
ensemble is reproduced with the choice: 



f(x) = e x exp 



F(x) 




x 2 ' 










x — a 


H 


- erf 









The energy constraint condition (|47|) in this ensemble is 



dyy^\A-y)e A y ' 



0. 



(64) 
(65) 

(66) 



The solutions A n ^ of this equation have been evaluated numerically for n = 
10, 20, 100, a = 2, 10, 20, and shown in Figured! For all the values of a investigated, 
A n ^ ~ |. The large n behaviour of A n ^ a is investigated in Appendix C The entropy, 
in the U representation is: 



S(U, V) = I log V + log Y n>a (U) + log c n J^-e% 



where 



Y n>a {U) = U$ / dxx% 



erf 



^4n,cr(l 



X 



a 



2a 



(67) 



(68) 
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n 



Figure 1. Solution of Eq. (|66|l . for different values of a and n 



In agreement with Eq. (1501 . the dependence of the entropy on U is of the type | \ogU . 
The integral Y n>a has been obtained by using the change of variable x = e/U, where e 
denotes energy 

The Gaussian ensemble of Eq. (1631) interpolates between canonical and 
microcanonical ensembles as does the Tsallis ensemble [11] . In fact, on one hand 
f(x) — > e x as cr — > +oo, and, on the other f{x) — > y/2TraS(x) as a — > (the vanishing 
term \/2ira is not a problem because it will be cancelled with the normalization). 

Further, based on the fact that the Gaussian ensemble of index a well describes 
the statistics of a large component of a large isolated system, we deduce that it 
must well approximate the Tsallis statistics of index a q = a in the case n, a q ^> 1. 
This fact is illustrated in Figure [2} where we have plotted the unnormalized Tsallis 

distribution ( 1 + 2^(1 _ f) J et _1 , for a q = 50, n = 60, U = 10, and fitted 

it to the unnormalized gaussian statistics exp[A(l — e/U) — i^-^/IDl jet -1 with 
U,o~,A as free parameters. The fit is very good (1 — R 2 = O(10~ 4 )) and the fitting 
parameter matched quite well the expected values: a = 48.36 ~ a q = 50, U = 10, 
A = 29.51 ~ A n=60>CT=50 = 29.63. 
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Gaussian fit of Tsallis statistics with a = 50, n = 60 and U = 1 

q 




Un-normalized Tsallis Statistics 
with a =50, n=60 
q 

- Fitted Gaussian 



FIT: 

y = exp[A(1-e/U)-A(1-e/U) 2 / (2a)] 
U = 10 a = 48.36 A = 29.51 



1-R 2 ~10" 4 



10 



15 



20 



25 



Figure 2. Gaussian fit of Tsallis statistics for a q — 50, n = 60, U = 10. The dual 
Gaussian ensemble approximates Tsallis ensembles for large values of n and a q 



5.3. Fermi-like ensemble 

The theory developed so far allows us to construct mechanical models of 
thermodynamics with the most diverse types of distributions. For example one may 
ask whether it would be possible to have an ensemble with a Fermi-like distribution. 
This is possible for the ideal gas case. We will describe this ensemble as an illustration of 
the theory, without discussing whether it really applies to some many-particle physical 
system. The Fermi-like statistics uses the choice: 



F(x) = hg(e x + 1) 
/(*) = l —r 

with no finite cut-off. The solutions A n of the energy constraint equation (J47J): 

y^- 1 (A-y) 



I(A,n) 



dy- 







(69) 
(70) 

(71) 



have been evaluated numerically and shown in Fig. [3] for the values n = 10, 20, 100. 



From the figure we see that, A n ~ n. This fact is analyzed in more details in | Appendix D 



Adopting the (U, V) parametrization, the ensemble, is: 



o(z: U. V) = 



+ 1 



n -i 



(72) 
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Figure 3. Solution of Eq. (fTTj) for different values of : 



or, alternatively: 



P(z;/W 



[ e -A n +f3H(z;V) + ij- 1 



(73) 



Jdz [e-A n +PH^y) + If 1 

if the (/?, V) parametrization is adopted. Using Eq. ( |26|) gives the entropy. Adopting 
the (U, V) parametrization, this is given by 



n 



S(U, V) = log c n + - log V + log J n (U) 



(74) 



where 



J n (U) = U* dxx^ 1 log (e A " (1 " x) + 1) (75) 
Jo 

In agreement with Eq. (150]) . the dependence of the entropy on U is of the type 
MogC/. The integral J n has been obtained by using the change of variable x = e/U, 
where e denotes energy. 



6. Discussion and conclusion 



In this work we addressed some fundamental issues raised recently in statistical 
mechanics, namely whether a theoretical basis can be provided for non-standard (i.e. 
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DUAL ENSEMBLES 




Others ? 



Micro canonical Canonical 



Figure 4. Generalization scheme. Tsallis ensembles, Gaussian ensembles, Fermi- 
like ensemble, and possibly other kinds of ensembles belong to the class of dual 
orthodes. Canonical and microcanonical are special instances of Tsallis and Gaussian 
ensembles. For large systems (n >> 1) Gaussian and Tsallis ensemble coincide if 
a = ^ = C v » 1. 

neither canonical nor microcanonical) ensembles, which are often encountered in the 
most diverse fields of physics. Following a line initiated in Refs. [21 [3], we used 
Boltzmann's original approach based on the "heat theorem", in order to provide a 
fresh look on the subject. By generalizing the duality property observed in the Tsallis 
case, we have been able to define the class of dual statistics, which includes the Tsallis 
ensembles as particular cases. The generalization scheme is represented in Fig. HI 
Thanks to the proposed generalization it is possible to provide a theoretical basis for 
many non-standard statistics other than Tsallis', such as the gaussian and the Fermi- 
like statistics. For all such non-standard orthodes the heat theorem holds, namely the 
usual thermodynamic relations are recovered. In this paper we also have provided a 
general formula for the entropy associated to any dual orthode. This allowed us to 
write down the expression of entropy for the Gaussian ensemble, which has never been 
done before. To use the same expression as Gallavotti [I], all dual orthodes provide 
"mechanical models of thermodynamics" . This result is not trivial since the only known 
orthodes were the microcanonical ensemble, the canonical ensemble and its variants like 
the grand canonical ensemble and the pressure ensemble [I]. Despite the fact that the 
class of orthodes is quite vast, the microcanonical and the canonical ensembles still play 
a special role in statistical mechanics because they are cases of "hidden dual statistics" . 
They do not rely on the employment of the energy constraint, which constitutes the 
mechanism through which it is possible to construct non-standard dual orthodes. 

One interesting result is that the classical ideal gas thermodynamics follows from 



Thermodynamics with generalized ensembles: The class of dual orthodes 



17 



any dual ensemble, thus revealing that this occurrence is not special for the standard 
statistics. This fact was already noticed for the Tsallis case a few years ago [10], although 
its orthodicity was not yet clearly recognized. It also suggests that there is a one-to-one 
correspondence of states obtained in different dual ensembles, i.e. there is an equivalence 
of all dual ensembles. Such equivalence holds no matter the number of degrees of 
freedom, and, of course, may break as one considers more realistic Hamiltonian models 
of systems with phase transitions. We also like to stress that all dual orthodes are 
equivalent also in another sense, namely the sense first investigated by Gibbs [12] for 
the canonical and microcanonical ensembles. In the thermodynamic limit (n — > oo) 
and for the free gas model, the canonical ensemble is so peaked around the average 
energy value that it is practically undistinguishable from the microcanonical one [13J. 
The same kind of equivalence occurs for any dual statistics provided that A n — > oo in 
the thermodynamic limit, which is indeed the case observed in the examples considered 
in this paper. This is because the distribution is expressed in terms of the quantity 
/ (A n ^p) which tends to the Dirac delta (times an unimportant proportionality factor), 
centered around U: 

f [An ~ ) -> const x 5(U - H) (76) 



U 

This follows from the asymptotic formula h(kx) — > dy ^ v ^ S(x), as k — » oo. 
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Appendix A. Proof of the Generalized Equipartition Theorem of Eq. (I32j) 

Let us calculate explicitly the average value of Pi^. Regardless of the parametrization 
one has: 



P 



1 



[ Pl F[(3(U - H(z; V))]]zem> - J dzF[P(U - H)] 



(A.l) 

Where we integrated by parts over pi to obtain the third line. From Eq.s ( 1271) and ( 1281 ) 
follows: 

/ 9H\ ig /A N 

According to the parametrization adopted this would be either Eq. ( 1301) or ( I3TT) . 
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Appendix B. Proof of orthodicity in the (3, V parametrization 



In the /3,V parametrization U is a function of (/3,V). Therefore F is a function of 
(z; U, V), and Q (Eq. [27]) and S p (Eq. [26]) are functions of (/?, V). Let us calculate the 
partial derivatives of S 



W = m 1 d * mu ~ H)) 

1 



dzf(P(U-H)) 



(U-H) + p 



1 dU n 



dU 



G ITJ __. GBU 

g {H ' u)f)+ g P dp = Y p ^p 



(B.l) 



where we used the equation {7 P =< > p from the state definition ([291) and the 
generalized equipartition theorem of Eq. ([32]) . 



as, 

<9K 



lJLJd*FWU-B)) 
dzf(P(U-H)) 



p gdv g p \ov/ i 
i 0f/„ p„ 



p dv p dv 



dH\ 



T p dV p 



T 



(B.2) 



where we used the expression for P p from Eq. (129]) and the generalized equipartition 
theorem of Eq. (1321) . Combining all together: 



(B.3) 



Appendix C. Large n behaviour of the coefficients A na 
Eq. (1661) can be recast in the following form: 



dyg n {y)e 2 - = A j dyg n (y)y e ^ 
o ./o 



where 



:*(!/) 



Jo ?/ 2e y 

and # is Heaviside's step function. Its Fourier transform is: 

Using the formula liniTv^oo (1 + 1 = e x , we have: 
lim £ n (Jb) = e ik ^ +1 ) 



(C.l) 

(C.2) 

(C.3) 
(C.4) 
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therefore: 

lim g n (y) = 8 (y-(^ + l)). (C.5) 
Using this result Eq. ( 1C.1I) . becomes, for very large n: 

A ~ | (C.6) 
This result does not depend on a. 

Appendix D. Large n behaviour of the coefficients A n for the Fermi-like case 

Figure [3] suggests that, in the limit n — ► oo, A n ~ n. In this appendix section we provide 
a simple consistency argument to support the claim that A n ~ n. Let us assume that 
for very large n, A n ~ n. Then we should have from Eq. (ITTT) : 

2/ 2 / , 2/ 2 



~ n / c/y (D.l) 

e^" n + 1 Jo e?/ ~ n + 1 

Equating the first derivative of the integrand in the left hand side to zero, gives: 

y = 2(1 + «,»-») (D.2) 

which is satisfied for y = n. This means that an extremum (a maximum as we will see) 

n 

is attained for y = n. The value taken by the integrand at the maximum is which 
increases very quickly. The second derivative, calculated at y — n, is: 
1 
4 

which tends to — oo very quickly. This indicates that the integrand becomes very sharply 
peaked around y = n as n increases. Therefore, as an approximation, we can replace y 
with n, and see that y% ~ ra/t -1 , thus getting Eq. (ID.lj) . 



71/ \ n_ _|_ i _ 

I n 2 — n 2 

2 



(D.3) 



References 

[1] C. Tsallis. Possible generalization of Boltzmann-Gibbs statistics. Journal of statistical physics, 
52:479, 1988. 

[2] Michele Campisi and Gokhan B. Bagci. Tsallis ensemble as an exact orthode. Physics Letters A, 

362(1):11-15, February 2007. 
[3] Michele Campisi. On the limiting cases of nonextensive thermostatistics. Physics Letters A, to 

appear, 2007. arXiv:cond-mat/0611068. 
[4] G. Gallavotti. Statistical mechanics. A short treatise. Springer Verlag, Berlin, 1995. 
[5] Ludwig Boltzmann. Uber die Eigenschaften monocyklischer und anderer damit verwandter 

Systeme. Crelle's Journal, 98:68-94, 1884. Reprinted in Hasenohrl (ed.), Wissenschaftliche 

Abhandlungen, vol. 3, pp. 122-152. New York: Chelsea. 
[6] M. Campisi. On the mechanical foundations of thermodynamics: The generalized Helmholtz 

theorem. Stud. Hist. Philos. M. P., 36:275-290, 2005. 
[7] A.I. Khinchin. Mathematical foundations of statistical mechanics. Dover, New York, 1949. 
[8] S. Martinez, F. Pcnnini, A. Plastino, and C. Tessone. On the equipartition and virial theorems. 

Physica A, 305, 2002. 

[9] M.P.Almeida. Generalized entropy from first principles. Physica A, 300:424-432, 2001. 



Thermodynamics with generalized ensembles: The class of dual orthodes 



20 



[10] S. Abe, S. Martinez, F. Pennini, and A. Plastino. Noncxtcnsive thermodynamic relations. Phys. 
Lett. A, 281:126, 2001. 

[11] Murty S. S.Challa and J. H.Hetherington. Gaussian ensemble as an interpolating ensemble. 

Physical Review Letters, 60(2):77-80, 1988. 
[12] J.W. Gibbs. Elementary principles in statistical mechanics. Yale University Press, Yale, 1902. 

reprinted by Dover, New York, 1960. 
[13] K. Huang. Statistical mechanics. John Wiley & Sons, Singapore, 2nd edition, 1963. 



